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ABSTRACT 

Spatial velocities of all currently known 28 masers having trigonometric parallaxes, 
proper motion and line-of-site velocities are reanalyzed using Bottlinger's equations. 
These masers are associated with 25 active star-forming regions and are located in 
the range of galactocentric distances 3 < i? < 14 kpc. To determine the Galactic 
rotation parameters, we used the first three Taylor expansion terms of angular rotation 
velocity Q. at the galactocentric distance of the Sun i?o = 8 kpc. We obtained the 
following solutions: fio = -31.0 ± 1.2 kms"^kpc"\ % = 4.46 ± 0.21 kms^^kpc^^, 
= -0.876±0.067kms"^kpc~^ Oort constants: A = 17.8±0.8 km s^^ kpc"\ B = 
— 13.2±1.5 km kpc~^ and circular velocity of the Solar neighborhood rotation Vb = 
248±14 kms~^. Fourier analysis of galactocentric radial velocities of masers Vr allowed 
us to estimate the wavelength A = 2.0 ± 0.2 kpc and peak velocity fn = 6.5 ± 2 km 
s~^ of periodic perturbations from the density wave and velocity of the perturbations 
4 ± 1 kms"""^ near the location of the Sun. Phase of the Sun in the density wave is 
estimated as x& ~ —130° ± 10°. Taking into account perturbations evoked by spiral 
density wave we obtained the following non-perturbed components of the peculiar 
Solar velocity with respect to the local standard of rest (LSR) {Uq,Vq,Wq)lsr = 
(5.5, 11, 8.5) ± (2.2, 1.7, 1.2) kms"\ 

Key words: Masers - SFRs - Spiral Arms: Rotation Curve - Galaxy (Milky Way). 



1 INTRODUCTION 

Study of Galactic parameters, especially of the Galactic ro- 
tation curve, is of great importance in solving a number of 
problems, such as estimating the mass of the Galaxy, de- 
termining the distribution of matter, estimating the hidden 
mass, studying dynamics and structure of the Galaxy and 
its subsystems, etc. 

For this purpose masers having trigonometric paral- 
laxes, line-of-sight velocities and proper motions (Reid et 
al., 2009) are of great interest. Currently measurements of 
parallaxes with mean error of several percents are fulfilled 
with radio interferometers for three tens of masers in the 
various active high-mass star-forming regions. Because these 
masers are associated with very young (< 10^ yr) OB stars, 
their kinematics reflects the properties of the youngest part 
of the Galactic disk. 

Analysis of motions of 18 masers was made by a number 
of authors (Reid et al., 2009; Baba et al., 2009; Bovy et al., 
2009; McMillan & Binney, 2010). For the first time Reid 
et al. (2009) suggested that high-mass star-forming regions 
significantly (~ 15 km s~^) lag circular rotation. McMillan 



& Binney (2010) suggest that much of this lag could be 
accounted for by increasing the value of the Solar motion in 
the direction of galactic rotation from 5 to 11 km s~^. 

The question of particular interest is the peculiar veloc- 
ity of the Sun with respect to the LSR, {Uq,Vq,Wq)-lsr. 
On the basis of Stromberg's relation, Dehnen & Bin- 
ney (1998) determined the vector {Uq,Vq,Wq)-ls'r = 
(10.0,5.3, 7.2)±(0.4, 0.6,0.4) kms"^ using proper motions of 
~ 12000 main sequence stars from the Hipparcos catalogue. 
These results were confirmed by Aumer & Binney (2009): 
(f/0,F0,lF0)LSR = (10.0,5.3,7.1) ± (0.3,0.5,0.3) kms-\ 
who applied the same method to proper motions from the 
revised version of the Hipparcos catalogue (van Leeuwen, 
2007). In the work of Schonrich et al. (2010), where gra- 
dient of metallicity of stars in the Galactic disk was taken 
into account, this velocity is different: {Uq,Vq,Wq)lsr = 
(11.1,12.2,7.3) ± (0.7,0.5,0.4) kms'V Using another 
approach, Francis & Anderson (2009) suggested that 
this velocity is: (f/g, T/q, 1F0)lsr = (7.5,13.5.2,6.8) ± 
(1.0,0.3,0.1) kms-\ 

In the present work, we are trying to establish relation- 
ship between motions of all currently known masers having 
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parallaxes, proper motions and line-of-sight velocities, and 
parameters of the Galactic spiral density waves (Lin & Shu, 
1964), and to estimate components of the peculiar velocity 
of the Sun with respect to the LSR. This goal is achieved 
by determining parameters of the Galactic rotation curve, as 
well as other kinematic parameters, by means of Bottlinger's 
equations. Fourier analysis of periodic deviations of observa- 
tional velocities of masers from the Galactic rotation curve 
is aimed to obtain some estimates of spiral density wave 
parameters. 



2 DATA 

Input data on methanol masers with trigonometric paral- 
laxes and proper motions from VLBI measurements are 
taken from the papers of Reid et al. (2009), Rygl et al. 
(2010). Wc added also H2O masers SVS 13 in NGC 1333 
(Hirota et al., 2008a), IRAS 22196+6336 in Lynds 1204G 
(Hirota et al., 2008b) and G14.33-0.64 (Sato et al., 2010), 
with parallaxes and proper motions measured in the frame- 
work of VERA (VLBI Exploration of Radio Astrometry) 
program. Note that, for mascr NGC 281 W listed in Reid 
et al. (2009), we used observational data obtained by Rygl 
et al. (2010) with the quoted uncertainty for 6.7 GHz masers 
smaller than the Sato ct al. (2008) 22 GHz masers. For the 
nearest star-forming region in Orion, along with the data 
on maser Orion-KL (Hirota ot al., 2007), wc use indepen- 
dent VLBA observations of radio star GMR A (Sandstrom 
et al., 2007) with line-of-sight velocity with respect to the 
LSR (K)lsr = -4 ± 5 km s"^ (Bower et al., 2003). It is 
the Bower's et al. (2003) opinion that the GMR A source 
is a young star of T Tau type (WTTS). We use also a very 
young object a supergiant S Per (Asaki et al., 2007) with 
line-of-sight velocity K- = —39.7 ± 0.7 kms^^ from Famaey 
et al. (2005). 

All 28 objects being studied are listed in Table 1. In 
Fig. 1, one can see only 25 points because X, Y coordinates 
of some objects are very close to each other. For example, 
three masers Cep A, IRAS 22198 and L 1206 form practically 
a single point in the XY plane. 

Fig. 1 is drawn in galactocentric Cartesian coordinates 
X and Y. As one can see from this picture, only one maser 
has inaccurately measured parallax. That is the maser W51 
(in the first Galactic quadrant), with relative error of par- 
allax en/n = 0.36. In fact, only this maser requires calcula- 
tion of upper and lower bounds on its distance. But, because 
this maser lies in the horizontal part of the Galactic rotation 
curve, hereafter we use an average estimate of the heliocen- 
tric distance r ± or galactocentric distance R ± cr_R. The 
same distance estimates are used also for other masers. 

We reduce line-of-sight velocities of masers with re- 
spect to the LSR to the heliocentric coordinate system 
using the standard Solar motion velocity: (a, 5)1900 = 
(270°,-F30°), V = 20 kma'\ hence {Ue,VQ,We)i^sn = 
(10.3, 15.3, 7.7) kms"^ (Reid et al., 2009). 

It is very important for the purposes of our study to 
know the accurate value of Ro- A review of Ro estimates 
made by a number of authors can be found in the paper 
by Avedisova (2005), where the average value is derived as 
Ro = 7.80 ± 0.33 kpc. Note that the most reliable values 
of this parameter were obtained only recently. Analysis of 




-1 1 2 3 4 5 
Y.kpc 



Figure 1. Coordinates of masers in the XY Galactic plane. Lo- 
cations of the Sun and Gala<;tic Centre (GC) are indicated. The 
circle of 8 kpc radius is drawn. 



motion of the star S2 around the black hole in the Galactic 
Centre gave Ro = 8.4±0.4 kpc (Chez et al., 2008). This value 
is in a good agreement with observations of orbits of 28 stars 
near the Galactic Centre during 16 yr: Ro = 8.31 ± 0.33 kpc 
(Gillessen et al., 2009). 

In general, modern estimations show that the most 
probable value of Ro is close to 8.0 kpc, and error of its 
determination is close to 0.3 kpc. For the sake of reliability, 
we'll consider that the value of Ro lies in the 7.5 — 8.5 kpc 
range. 



3 METHODS 

3.1 Determining Rotation Curve Pcirameters 

Method used here is based on the well-known Bottlinger's 

formulas (Ogorodnikov, 1958), where angular velocity of 
Galactic rotation is expanded in a series to n-th order terms 
in r/Ro: 

Vr = —uq cos b COS I — vq cos 6 sin i — wq sin b — (1) 
-Ro sin I cos b[{R - Ro)nl/V. + ... + {R- RorQ^/n\], 
Vi = uq sin I — vq cos I + rOo cos b — (2) 
-{Ro cos/ -rcos b)[{R- Ro)il'o/l\+ 
+ ... + {R-R^rfl^/n\\, 

Vb = Uq cos I sin b + VQ sin Z sin 6 — wq cos b + (3) 
+Ro sin I sin b[{R - _Ro)Oj/l! + ... + {R- RoTnl^ /n\] , 

where Vr is heliocentric line-of-sight velocity; V; = 
4.74T7i(Cos6 and Vt, = 4.74r/Xj, are proper motion velocity 
components in the / and b directions, respectively (the co- 
efficient 4.74 is the quotient of the number of kilometers 
in astronomical unit by the number of seconds in a tropi- 
cal year); r = I/tt is the heliocentric distance of an object; 
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. Data, on inasors. 


Volocitios (/^. 


.11 ai'(~ rt^sif 


ual volocitios 


witli Galactic rotation oxcliidcd using parairiotoi> 


(16 


Source 


R, kpc 


U^, km 


yl, km s-l 


W'^, km s-i 


AVe, km s-i 


Vr, km s 1 


Ref 


L 1287 


8.52 


—2.9 di 1.7 


on T _L o a 
— 20, ( ± 2.D 


— 10.2 ± 2.5 


—5.7 ± 2.6 


—5.6 ± 1.7 




IKAb UU42U+5ooU 


9.33 


3.0 lb 2.7 


— 19.6 lb 4.2 


— 1.5 ± 0.8 


—2.9 ± 4.2 


-11.8 ± 2.8 


(1) 


JNGC 2ol— W 


9.50 


—4.2 lb 2.6 


—2.3 lb 2.8 


— 16.0 ± 1.7 


12.8 ± 2.8 


— 1.2 ± 2.6 


(2) 


W3 (OH) 


9.46 


4.8 lb 2.1 


— 17.D ± 2.2 


—6.2 ± 0.2 


— 1.2 ± 2.2 


— 13.1 ± 2.1 


(1) 


WB 89-437 


13.05 


0.2 lb 2.1 


— 17.1 ±2.1 


2.5 ± 3.8 


0.2 ±2.1 


—8.6 ± 2.1 


(1) 


NGO 1333 (rlj 


8.21 


— 23.D ± 4.4 


— 13.5 ± 2.3 


0.0 ± 2.5 


0.9 ± 2.3 


15.6 ± 4.4 


(3) 


JNGO 1333 (i2) 


8.21 


—21.1 ± 4.5 


5.5 ± 1.8 


10.0 ± 3.5 


20.0 ± 1.8 


13.3 ± 4.5 


(3) 


Un CjrMH A 


8.32 


—4.7 ±4.1 


—8.8 ± 2.3 


—3.1 ± 1.7 


5.6 ± 2.3 


—3.4 ± 4.1 


(4) 


Ori KL 


8.34 


— 17.9 it 4.2 


-\ A A 1 rt A 

— 14.4 ± 2.4 


—3.8 ± 3.0 


0.3 ± 2.4 


9.9 ± 4.2 


(1) 


S252 A 


10.08 


— 12.D ± 3.0 


—11.1 ± 0.6 


—9.2 ± 0.3 


3.5 ± 0.6 


4.5 ± 3.0 


(1) 


S255 


9.56 


—7.7 ± 3.6 


—0.8 ± 9.5 


—4.0 ± 10 


13.7 ± 9.5 


—0.8 ± 3.6 


(2) 


S269 


13.16 


—3.5 ± 2.9 


— 16.2 ± 1.2 


-11.7 ± 0.9 


—2.2 ± 1.2 


—4.3 ± 2.9 


(1) 


V Y OMa 


8.64 


—6.1 ± 2.9 


— 16.9 ± 3.0 


— 13.6 ± 2.5 


—2.7 ± 3.0 


— 1.6 ± 2.9 


(1) 


Gi32. 62+0.9 


9.12 


—7.8 ± 3.2 


— 10.5 ± 3.1 


—6.4 ± 2.2 


4.0 ±3.1 


—0.9 ± 3.2 


(1) 


(jl4. 33— 0.64 


6.92 


it 10 


— 10 ± 14 


— 10.8 ± 4.5 


5 ± 14 


—8 ± 10 


(1) 


Cj23. 06—0.13 


5.23 


21. o it 3.7 


—6.7 ± 5.7 


Q 1 -i- n yl 
— o.l It U.4 


14. o It O.D 


OK Q -I- Q C 
— ZD.O It O.O 


(1) 


(jr23.Ul — U.41 


4.18 


1.0 it 4.7 


± 9.1 


-8.5 ±2.3 


-1.8 ±8.4 


-10.8 ±5.8 


/I ^ 

(1) 


VjrZO.'*0 — U. J-O 


o.ou 


1 -1- 1 n 

— ±o in J-U 


Q -1- 99 


-5.3 ± 1.8 


13 ±18 


20 ±16 


(^^ 


G35.20-0.74 


6.34 


-21.7 it 3.3 


-20.1 ±3.6 


-15.7± 1.5 


-8.2 ±3.6 


12.2 ±3.3 


(1) 


W48 


5.65 


-31.1 it 5.7 


-16.5 ±7.5 


-16.5 ±2.4 


-9.6 ±7.3 


21.1 ± 5.9 


(1) 


W51 


6.08 


-24 ±41 


-10 ±35 


-10.4 ±3.7 


-7 ±37 


15 ±38 


(1) 


V645 


7.16 


-14.4 ±2.6 


-13.6 ± 2.9 


-11.4 ±0.5 


-0.8 ± 2.9 


6.3 ± 2.6 


(1) 


Onsalal 


7.50 


-9 ± 13 


-18.3 ± 5.6 


-2.9 ± 6.0 


-4.1 ± 6.8 


± 12 


(2) 


IRAS 22198 


8.26 


-4.2 ±3.7 


-23.4 ±4.9 


3.0 ±2.1 


-8.6 ±4.9 


-4.6 ±3.7 


(3) 


L 1206 


8.27 


-12.8 ±4.4 


-19.4 ±3.2 


-7.2 ±5.7 


-5.3 ±3.2 


4.4 ±4.3 


(2) 


Cep A 


8.26 


-8.9 ±3.8 


-16.6 ±4.9 


-12.6 ± 1.2 


-2.2 ±4.8 


0.7 ±3.8 


(1) 


NGC 7538 


9.30 


-2.1 ±2.1 


-32.7 ±2.9 


-18.1 ±1.0 


-16.0 ±2.8 


-10.5 ±2.2 


(1) 


S Per 


9.30 


-16.1 ±1.2 


-16.9 ± 1.2 


-17.2 ±3.1 


-3.9 ± 1.1 


7.5 ±1.1 


(5) 



Notes. Initial data are taken from: (1) - Reid et al. (2009); (2) - Rygl ct al. (2010); (3) - Hirota et al. (2008a); (4) - Sandstrom et al. 

(2007); (5) - Asaki et al. (2007). 



proper motion components /ij cos 6 and iih are in masyr"^, 
line-of-sight velocity Vr is in kms~^; «©,«©, w© are Solar 
velocity components with respect to the mean group veloc- 
ity under consideration; _Ro is the galactoccntric distance of 
the Sun; R is the galactocentric distance of an object; 7?, 
Rq and r are in kpc. Velocity U directed towards the Galac- 
tic Centre, V along the Galactic rotation, and W towards 
the Northern Galactic pole. The quantity Qo is the Galac- 
tic angular rotational velocity at distance Ro, parameters 
, • • • 5 ^^0 s^re derivatives of the angular velocity from the 
first to the n-th order, respectively. The distance R can be 
calculated using the expression 



R' 



cos^ 6 — 2 J?or cos h cos i + i?o • 



(4) 



The system of conditional equations (l)-(3) contains n + 4 
unknowns: Uq,vq,Wq, f2o, ^J, . . . , f2o , which can be deter- 
mined by the least-squares method. The system of equations 
(l)-(3) is solved with weights of the form 



(r,l,b) 



(5) 



where Pr,Pi and Pb are the weights of equation for Vr, Vj 
and Vb correspondingly, and So = 8 kms~^ is the "cosmic" 

dispersion averaged over all observations. Errors avi and av^ 
of velocities Vi and Vb can be calculated using the formulas 



(6) 



3.2 Determining Projections Ve and Vr 

Components of spatial velocities U, V, W of masers are de- 
termined from observed radial (line-of-sight) and tangential 
velocities Vr, Vi, Vb in the following way: 

U = Vr cos I cos b —Vi sin I — Vb cos / sin 6, 
V = Vr sin I cos b + Vi cos I — Vb sin / sin b, (7) 
W — Vr sin b + Vb cos b. 

Then we find two projections of these velocities: Vr, directed 
radially from the Galactic Centre towards an object, and Ve, 
orthogonal to Vr and directed towards Galactic rotation: 



Ve = Usme + {Vo + V) cose, 
VR = -Ucose + (Vo + V) sin 6, 



(8) 
(9) 



where Vo = |i?()0()|, and position angle 9 is determined as 
tan6' — y/(J?() —x), where x, J/ are Galactic Cartesian coordi- 
nates of an object. In addition, it is assumed that velocities 
U and V are free from the Solar velocity with respect to 
mean group velocity (uq,vq,wq) obtained from equations 
(l)-(3). Errors of projections Vr and Ve are estimated from 
the expressions 



2/T7- \ 2 ■ 2 a I 2 2 , 

a [Ve) = au sm. + ay cos I 

2/TT \ 2 2 /] I 2.2 

f" (Vr) = COS 6^ + <T^ sm 



(10) 

(11) 

where errors of spatial velocities U and V are denoted by 
(7(7 and av, respectively. 

We are interested also in heliocentric spatial velocities 
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7 8 9 10 11 12 13 14 
R, kpc 

Figure 2. Galactic rotation curves Vg{R) for adopted Ro 
7.5,8.0,8.5 kpc, lines 1,2 and 3. 



which are free from Galactic rotation with parameters (16) 
but not from Solar velocity. We denote such velocities by 
[/\ y\ (sec Table 1 and Fig. 7). These velocities must 
satisfy the following equations: 

u' = -ue, v' = -ve, W = -wq. (12) 



3.3 Fourier Analysis of Velocities 

At the next step, we consider deviations from circular veloc- 
ities AVe and galactocentric radial velocities Vr. 

Spectral analysis consists in applying direct Fourier 
transform to the sequence of residual velocities following 
way: 



JT-Ri 



(13) 



where V{Xk) is the fe-th harmonic of Fourier transform, M 

is the number of measurements of velocities with coordi- 
nates Ri, 4 = 1,2,..., M, and Xk is wavelength in kpc; the 
latter is equal to D/k, where D is the period of the origi- 
nal sequence in kpc. In density wave theory, A denotes the 
distance between adjacent spiral arms along the Galactic 
radius vector. 

Understanding the results of spectral analysis is easy in 
the framework of linear theory of density waves. According 
to this approach (Burton & Bania, 1974; Byl & Ovenden, 
1978), 



Vr = -/ijCOSX, 
AVe = /esinx, 



(14) 
(15) 



where x = 'rn[cot{i) \n{R/Ro) — 0]+Xe is phase of the spiral 
wave (m is number of spiral arms, i is pitch angle, xo is the 
phase of the Sun in the spiral wave (Rohlfs, 1977)); Jr and 
fe are amplitudes of radial and tangential components of 
the perturbed velocities which, for convenience, are always 
considered positive. 

For the case when perturbations Vr and AVe are 
considered as periodic functions of logarithm of R, the 
following expression for the phase is appropriate: x = 
(27ri?o/A) ln{R/Ro) + x@ ^ (27r/A)(i? - Ro) + Xo (Mel'nik 
et al, 2001), what allows us to use Fourier analysis (13). 

In our coordinate system radial velocity Vr is directed 




2 3 4 5 6 7 8 9 10 11 12 13 14 
R, kpc 

Figure 3. (a) Rotational angular velocity and (b) rotational cir- 
cular linear velocity of the Galaxy according to (16) versus galac- 
tocentric distance. Dashed lines denote the bounds of the confi- 
dence interval corresponding to the la level. Vertical line indicates 
the location of Ro = 8 kpc. 



from galactic centre and circular velocity deviation AVe di- 
rected along Galactic rotation. Centre of the spiral arm cor- 
responds to phase x = 0. For dcfiniteness (Rohlfs, 1977), 
we consider that the centre of the Carina-Sagittarius arm 
{R !^ 7 kpc) corresponds to phase X = 0: while the Perseus 
arm {R m 9 kpc) corresponds to phase x ~|^- 



4 RESULTS 

4.1 Galactic Rotation Curve 

Firstly we investigated solutions of equations (l)-(3) using 
the different number of Taylor expansion terms of fl for 
different values of Rq. Direct construction of the rotation 
curve using the rotation parameters obtained above showed 
that every additional expansion term leads to significant 
widening of the confidence interval. We revealed that it is 
enough to take three terms of Taylor expansion to construct 
a smooth, sufficiently accurate Galactic rotation curve in 
the galactocentric distance range of 3 < i? < 14 kpc. In- 
deed, adding subsequent terms to the expansion practically 
does not change the shape of the curve, and additional terms 
do not differ statistically from zero. So, in this sense we 
found the optimal solution. Results of this study are shown 
in Fig. 2, whence wc can see that the shape of the rota- 
tion curve significantly depends on the adopted value of 7?o. 
Thus, variation dRo = 0.5 kpc leads to variation of the ro- 
tational angular velocity dVe w 20 kms~^. As one can see 
from Fig. 2, influence of inaccuracy in Ro is particularly 
strong in the Solar neighborhood. 

Then we obtained solution of (l)-(3) for the fixed 
value of i?o = 8.0 kpc for 84 equations. As a result we 
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Figure 4. Deviations from circular velocities AVe after subtrac- 
tion of the model rotation curve with parameters (16) versus 
galactocentric distance; vertical line indicates Rq = ?> kpc. 



found the following parameters of the Solar velocity with re- 
spect to the group velocity: {uq, vq, wq) = (8.0, 14.9, 7.6) ± 
(2.1, 1.9, 1.6) kms"'^, and parameters of the angular velocity 
of Galactic rotation: 



Qo = -31.0 ± 1.2 kms"^kpc"\ 

^0 = +4.46 ± 0.21 kms"^kpc"^ 
nf,' = -0.876 ± 0.067 kms"^kpc"^. 



(16) 



The unit weight error is ctq = 9 km s . Oort constants A = 
0.5^0^0, B = no+0.5Jiofio are: A = 17.8±0.8 kms'^kpc"^ 
and B = —13.2 ± 1.5 kms~^kpc~^. Using these values 
and taking into account the error aRg = 0.3 kpc, we es- 
timated the circular velocity of the Solar neighborhood 
Vb = |-Rof2o| = 248±14 kms~^ and its orbital period around 
the Galactic Centre T = 27r/fio = 198 ± 11 Myr. 

Galactic rotational angular velocity Q{R) and the cor- 
responding linear circular velocity Ve{R) ~ \R^l{R)\ curves 
are shown in Figs. 3a and 3b, respectively. 

As one can see from Fig. 3a, errors in parameters (16) 
lead to strongly increasing errors in the rotational angular 
velocity Q{R) when _R > 10 kpc. Due to this the confidence 
interval for the linear rotation curve Vg{R) (Fig. 3b) starts 
to exceed the criterion adopted above, dVg = 20 km s^^, at 
ii > 10 kpc. In general, we can conclude that the Galactic 
rotation curve found above approximates our initial data 
with fully adequate accuracy over the range of distances of 
3 < i? < 10 kpc. 



4.2 Periodic Perturbations 

In Fig. 4, we present distribution of deviations from circu- 
lar velocities AVg versus R obtained by subtraction of the 
rotation curve found above, with parameters (16), from the 
initial data. Fourier power spectra estimates (periodograms) 
of the sequence AVe (|AVe|^) and galactocentric radial ve- 
locities Vr (|T^|2) versus wavelength A are shown in Figs. 5a 
and 5b, respectively. In case of AVg, using two main low- 
frequency Fourier spectrum peaks extending up to A = 

1.3 kpc, we constructed an approximation curve of the resid- 
uals, as shown in Fig. 4. The amplitude of the approximation 
curve was fitted using a chi-square optimization technique. 
As it is seen from the picture, the approximation curve is 
not representative enough to make any conclusions. 

But we have completely different situation in the case 
of radial velocities versus -R shown in Fig. 6. In accordance 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 
X. kpc 

Figure 5. Periodograms of (a) deviations from circular velocities 
AVg and of (b) radial velocities Vr. 



with the main periodogram peak shown in Fig. 5b, which 
is the largest one in the frequency region of interest, per- 
turbations in radial velocities have an average wavelength 
of A = 2.0 ~ 0.25i?o kpc, which is in good agreement with 
results of other authors (see section 5.2). The approxima- 
tion curve built using this main low-frcqucncy Fourier spec- 
trum peak extending from A = 3 kpc to A = 1.5 kpc and 
fitted using the chi-square optimization, is shown in bold 
in Fig. 6. One may clearly see that this curve gives the 
most exact description of radial velocities in the region of 
the Perseus arm (i? ~ 9 kpc). In this region, extended as 
far as the location of the Sun, errors of amplitudes of ra- 
dial velocities are less than ±1 kms~^. Peak value of the 
approximation curve reaches Jr = 6.5 kms"'^ and at the 
area close to the Sun it is about 4.2 km s~^. Confidence in- 
tervals for the approximation curve were determined using 
Monte-Carlo simulation for 1000 random samples. In agree- 
ment with the approximation curve, phase of the Sun lies in 
the interval —140° < x© < —120° with the mean value of 
X© ~ —130°. Radial perturbations in the Solar vicinity we 
estimate as periodic function with wavelength A = 2 ± 0.2 
kpc and amplitude /r = 6.5 ± 2 kms~^. 

As it is seen from Fig. 5a and Fig. 5b both periodograms 
reveal peaks at wavelength A = 2 kpc. Wo estimated prob- 
abilities p of belonging of these peaks to the periodic sig- 
nal (not to noise) using simple formula derived by Vityazev 
(2001): |F(A)|Vc7g/M ^ -ln(l -p), where |F(A)p is pe- 
riodogram peak value at wavelength of interest, ctq is dis- 
persion of input data sequence. Ad is number of data. In 
the first case we obtained p = 0.6 and in the second case 
p = 0.992 what justified more reliable estimation of spiral 
wave parameters from Vr rather than from AVe- 



6 V. V. Bobylev & A.T. Bajkova 



V 10 





- j. 













2 3 4 5 6 7 8 9 10 11 12 13 14 
R, kpc 

Figure 6. Radial velocities of masers Vr versus galactocentric 
distance R. Thin curves mark the confidence interval correspond- 
ing to 1(7 level. Vertical line indicates the galactocentric distance 
of the Sun Ro = 8 kpc. 

4.3 Peculiar Solar Velocity with Respect to the 
LSR 

In Fig. 7, spatial velocities and (see Ta- 

ble 1), free from Galactic rotation according to pa- 
rameters (16), are given. Residual velocity compo- 
nents averaged over 28 masers are u ,v ,w = 
(-8.6,-13.6,-7.5) ± (2.1, 1.6, 1.3) kms~\ with dispersion 
vector (aui,aYi,cT^7i) = (10.8,8.7,6.8) kms~^. 

As one may see from this figure, residual velocities of 
most of masers form a sufficiently compact ellipse of smaller 
dispersion. Most of deviations from this compact distribu- 
tion are due to the following three masers: G23.66— 0.13, 
G23.43-0.18 and NGC 1333 (f2). The detail NGC 1333 (12) 
reflects inner velocity dispersion within maser NGC 1333 
which itself forms an expanding shell (Hirota et al., 2008a). 
Maser G23.66— 0.13 has significantly negative radial velocity 
which reflects a certain peculiarity in its motion. 

Excluding these three sources, for 25 masers forming 
a compact ellipse, we deduced the following mean veloc- 
ity [7\F\Ty^ = (-9.1,-15.5, -8.5)±(1.9, 1.3, 1.2) kms"\ 
with dispersion (o-^i , cr^i , o■^^-l ) = (9.3,6.6,6.1) kms~^. 

The question of great interest is how density wave per- 
turbations affect the values of the LSR velocity components. 
Let us denote components of "non-perturbed" peculiar Solar 
velocity with respect to the LSR as {Uq,Vq^Wq)-ls'r- For 
the local case, i.e. for the immediate Solar neighborhood, 
under the condition x ~ X© have 

U' = -(f70)LSR - Vb cos e + AVe sin 9, (17) 

= -(V'0)lsr + Vr sine + /We cos 6). (18) 

When 9 ^ Q {R ^ Ro), we have sinS — >• and cos 61 — >• 1. 
Then 

= -(?7©)LSR + /flCosx©, (19) 

= -(V0)lsr + fe sinx©, (20) 

where fR,fe > 0. As one can see from (19)-(20), signs 
of corrections caused by spiral density waves are deter- 
mined only by signs of cosx© and sinx©- Adopting that: 
mean amplitudes of spiral density wave perturbations in 
both tangential and radial velocities are 6.5 ± 2 kms^^ and 
X© = -130° ± 10°, we obtain, according to (19)-(20), the 



following "non-perturbed" velocity of ((7q, V©, H^q)lsr = 
(5.5, 11,8.5) ± (2.2, 1.7, 1.2) kms"\ This result agrees with 
the estimate of Solar motion by Schonrich et al. (2010) and 
Francis Sz Anderson (2009). 

We could not determine the value of x© using deviations 
from circular velocities of masers. But in order to confirm our 
results obtained above using Galactic radial velocities Vr of 
masers, we analyzed also a sample of Open Star Clusters 
(OSCs) from the compilation of Bobylev et al. (2007) which 
was based on the catalog of Piskunov et al. (2006), along 
with the line-of-sight velocities from the CRVOCA catalog 
(Kharchenko et al., 2007). We selected 117 young OSCs (age 
< 50 Myr) in radius r < 2 kpc using the following restriction 
for modulus of residual velocities: V AU^ + AV'^ + AW'^ < 
50 kms~\ 

Results of analysis of deviations from circular velocities 
AVg using Galactic rotation parameters (16) and galacto- 
centric radial velocities Vr are shown in Figs. 8a and 8b, 
respectively. We can see a good agreement between the 
approximation curves shown in panels (a) and (b): wave- 
lengths are almost similar, difference in phase is about 
7r/2, peaks are also comparable (6-8 km s~^). Periodograms 
of deviations from circular velocities and radial velocities 
are shown in Fig. 9. To obtain approximation curves in 
Fig. 8, we used low-frequency Fourier spectra up to wave- 
lengths A = 1.3 kpc and fitted the amplitude of pertur- 
bations using the chi-square optimization. Spectral anal- 
ysis of Vr (Fig. 9b) reveals periodic motion with wave- 
length A = 1.9(0.23i?o) kpc. Spectral analysis of AVg 
(Fig. 9a) displays three dominant peaks at wavelengths 
A = 2.21(0.28i?o),l.l(0.137?o),0.47(0.058i?()) kpc (compare 
with results by Clemens (1985):A = 0.22i?o, 0.12i?o, 0.05iio). 

Now we can compare these results with those obtained 
using radial velocities Vr of masers. We can see a good agree- 
ment in the behavior of the approximation curves shown 
in Fig. 8b and Fig. 6 in the Solar area. Fourier spectrum 
of radial velocities of masers Vr reveals a dominant peak 
of amplitude 6.5 kms~^ at wavelength 2 kpc. In addition, 
there is a good agreement of phase x© ^ —120° using 
data shown in Fig. 8b with the result obtained from masers 
X© « -130° ± 10°. 



5 DISCUSSION 

5.1 Galactic Rotation Curve 

The shape of the Galactic rotation curve obtained above 
is in a good agreement with the one obtained by Zabolot- 
skikh et al. (2002) for Ro = 7.5 kpc as a result of thor- 
ough reconciliation of distance scales of different star groups. 
The value f2o ~ —31.0 ± 1.2 kms^^kpc^^ (see solution 
(16)) is in a good agreement with the following result 
obtained from 18 masers: -29.9 kms"^kpc"^ < flo < 
-31.6 kms"^kpc"^ (McMillan & Binney, 2010) and with 
Qo = -30.3±0.9 kms^^kpc"^ (Reid et al., 2009). Our value 
Vb = 248 ± 14 kms~^ {Ro = 8 kpc) is in a good agreement 
with 254 ± 16 kms"^ {Ro = 8.4 kpc) by Reid et al. (2009) 
and 244 ± 13 kms"^ {Ro = 8.2 kpc) by Bovy et al. (2009). 
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Figure 7. Heliocentric spatial velocities and free from 
Galactic rotation with parameters (16). 



5.2 Spiral Structure Parameters 

The estimate A = 0.22_Ro of the distance between ad- 
jacent spiral arms was obtained by Clemens (1985). For 
Ro = 8.5 kpc used here it gives A = 1.8 kpc which is in 
a good agreement with the value A = 2.0 kpc obtained in 
our study using 28 selected masers, as well as with A = 1.9 
obtained from 117 young OSCs. Mel'nik et al. (2001) ob- 
tained the estimate A = 2.0 ±0.2 kpc by analyzing velocities 
of OB associations in the 3-kpc Solar neighborhood. An in- 
dependent estimate A = 1.7 ± 0.5 kpc was obtained from 
the analysis of Vr velocities of young OSCs in the distance 
range of 6 < < 9 kpc (Bobylev et al., 2008). 

Prom the analysis of OB stars and Cepheids, 
Byl & Ovenden (1978) found the following parameters of 
the Galactic spiral structure /_r = 3.6 ± 0.4 kms~^ and 
}e = 4.7 ± 0.6 kms-\ xo « 165 ± 1°, i = -4.2 ± 0.2°. 
According to Clemens (1985), the value of perturbation am- 
plitude is \ fe\ ~ 5 kms~^ at the wavelength of A = 0.22_Ro. 

Using Cepheids, Mislmrov et al. (1997) found the fol- 
lowing parameters: fa = 6.3 ± 2.4 kms^^ and fe = 4.4 ± 
2.4 kms-\ xo « 290 ± 16°, i = -6.8 ± 0.7°. Using Hippar- 
cos OB stars and Cepheids, Fernandez et al. (2001) showed 
that amplitudes fu and fg are about 6 — 4 kms~^ and 
20° < X© < 284°. Mean amplitude of radial perturbations 
/h = 6.5 ± 2 kms~^ found by us using masers is in a good 
agreement with the estimates listed above. 

Location of the Sun in the external part of the Carina- 
Sagittarius arm is in agreement both with the distribution of 
stars and gas (Russeil, 2003) and with density wave theory. 



5.3 Peculiar Velocity of the Sun 

It is worth mentiouiug some results of determination of the 
non-perturbed LSR velocity parameters obtained simultane- 
ously with parameters of the Galactic rotation, taking into 
account perturbations evoked by the spiral wave density. 
For example, these are: {Uq, Vq)lsr ~ (2, 5) ± (2, 2) kms^"'^ 
(Mishurov et al., 1997); (?70,Vq)lsr = (7.8,13.6) ± 
(1.3, 1.4) kms^i (Mishurov & Zenina, 1999); {Uq,Vq)i,sr = 
(9,12) ± (1,1) kms"^ (Lepine et al., 2001). Let us men- 
tion also the results of Fernandes et al. (2001), obtained 
using different samples of OB stars: (f/©, V©, H^0)lsr = 
(8.8,12.4,8.4) ± (0.7,1.0,0.5) kms"^ and of Cepheids: 
{Uq, Vq, W^0)lsr = (6.5, 10.4, 5.7) ± (1.2, 1.9, 0.7) kms'^. 
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Figure 8. Deviations from (a) circular velocities AVg and (b) ra- 
dial velocities Vr for 117 young OSCs versus galactocentric dis- 
tance R. Approximation curves axe shown with solid lines. Verti- 
cal line indicates galactocentric distance of the Sun Rq = 8 kpc. 
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Figure 9. Periodograms of deviations from (a) circular velocities 
AVg and (b) radial velocities Vr for 117 young OSCs. 



From the results listed here, we can see that values of 
([/, V) components of the LSR velocity are always smaller 
than the corresponding components of the standard Solar 
motion velocity {Uq, Vq, Wq) = (10.3, 15.3, 7.7) kms"^ 



6 CONCLUSIONS 

Spatial velocities of all currently known 28 masers hav- 
ing estimates of trigonometric parallaxes are reanalyzed. 
To determine the Galactic rotation curve, we used the 
first three terms of Taylor expansion of the Galactic an- 
gular rotational velocity at the galactocentric distance of 
the Sun Rq = 8.0 kpc. Coefficients obtained from Bot- 
tlinger's equations are: fio = —31. Oil. 2 kms~^kpc~^, fig = 
4.46±0.21 kms-^kpc"^, Q,'o = -0.876±0.067 kms'^kpc"^. 
Oort constants are A = 17.8 ± 0.8 kms~^kpc~^ and B = 
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— 13.2 ± 1.5 kms ^kpc ^, which allowed us to estimate 
the circular velocity of the Solar neighborhood rotation 
Vb = 248 ± 14 kms"^ {Ro = 8 kpc). 

It was found that deviations from circular velocities 
AVg of masors arc not representative enough for any reliable 
analysis, while the study of galactocentric radial velocities 
Vr of masers allowed us to estimate a number of useful pa- 
rameters concerning spiral structure of the Galaxy. Based on 
Fourier analysis of velocities Vr, we obtained that: there are 
periodic perturbations evoked by the spiral density wave, 
with wavelength of A = 2.0 ± 0.2 kpc and peak velocity of 
/r = 6.5 ± 2 kms~^ and velocity of 4 ± 1 kms~^ near the 
location of the Sun; phase of the Sun in the density wave is 
Xq = —130° ± 10°, which proves that the Sun is located in 
the inter-arm space. 

Adopting the parameters of spiral density wave found 
we determined the following non-perturbed Solar pecu- 
liar velocity with respect to the LSR: {Uq,Vq,Wq)lsr = 
(5.5,11,8.5) ± (2.2,1.7,1.2) kms"\ In this connection, in 
the regions near the Sun, velocity component (V0)lsr is 
influenced mostly by negative perturbation AVg and has 
smaller value than the corresponding velocity component 
Vq = 15.3 kms~^ of the Standard Solar motion. 
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